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A search for CP(T) violation is performed and the fractional difference between the b and b 
hadron lifetimes is measured using reconstructed secondary vertices in inclusive B hadron decays 
selected from 3.1 million Z° — > qq events. The data were collected by the OPAL experiment at the 
LEP collider at CERN at ~ 91 GeV from 1991-1995. 



I. INTRODUCTION 



CP violation in the B-meson system has generated considerable experimental and theoretical interest, as potentially 
large effects are expected. Searches for CP(T) violation using the small sample of Z° — > bb decays at the LEP collider 
at CERN provide "proofs of principle" for analysis techniques which will be employed by future B-Factory experiments. 

Indirect CP violation is possible in the B° system provided the weak eigenstates B° and B° differ from the mass 
eigenstates Bi and B 2 : 



|Bi > = 



(1 + e B + 5 B )\B° > + (1 - e B - 8 B )\B° > 



|B 2 > 



1 + \e B + S B 
(1 + e B - S B )\B° >-(l-e B + S B )\B° > 



2(l + \e B -S I 



(1) 



(2) 



where e B and S B parametrize CP and CPT violation, respectively J]J]. These parameters have been investigated using 
semileptonic b hadron decays, resulting in limits of order 10~ 2 on both e B and S B |^J^,^,||. In the Standard Model, 
Re(eB) is expected to be around 10~ 3 (6), but it could be up to an order of magnitude larger in superweak models 0. 

A non-zero value of es gives rise to a time-dependent rate asymmetry, A(t), in inclusive B° vs. inclusive B° decays 
H, defined as: 



A(t) 



B°(t) -B°(t) 
B°(t)+B°(t) : 



(3) 



where B°(t) and B°(t) are the decay rates of B° and B° mesons. For an unbiased selection of B° and B° mesons, the 
time-dependent inclusive decay rate asymmetry can be rewritten in terms of proper decay time t: 



A{t) = a c 



Am d T B a . 2 ( Am d t 
sm(/\mdt) — sin 



(4) 



where a cp is the CP-violating observable, Amj is the B° oscillation frequency, and tbo is the B° lifetime. 
For \e B \ << 1, the parameter a cp is related to e B by: Re(eB) = a cp /4. 

Furthermore, CPT invariance implies that Tf, = tj. If CPT violation occurred, the lifetimes of b and b hadrons 
could be different: 



T b/b 



1 

1 ± - 

2 



At 
r 



(5) 



where r av is the average and (Ar/r)b is the fractional difference in lifetimes. 
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II. INCLUSIVE CP(T) TESTS 



The measurement of the time-dependent rate asymmetry, A(t), and the extraction of Re(es) proceeds in several 
steps. First, selected Z° — > qq events are divided into 2 hemispheres defined by the plane _L to thrust axis and 
containing the e + e~ interaction point. A sample of about 400,000 Z° — > bb events is identified using b-tagging 
techniques described in detail in |^|,p| JIol[Tl|| . The b-tags rely on the presence of a displaced secondary vertex or a high 
momentum lepton. In each event, the hemisphere containing the b-tag is referred to as the "T-tagged" hemisphere. 
The b-tag has a ehcmi of 37 % for bb events and a non-b impurity of 13%. Next, the b hadron proper decay time, 
t, in the opposite "measurement" hemisphere is reconstructed by forming a secondary vertex, measuring the decay 
distance from the primary vertex, and estimating the b hadron energy. The quantity a cp is then extracted via a 
binned x 2 -fit to the observed time-dependent asymmetry in bins of reconstructed proper time. 



A. Production Flavor Tag 

The production flavor estimate, Qt, in the tagged hemisphere is the output of a neural net with the following 
inputs: 

1. Jet charge of the highest energy jet, Qj C t, with (with momentum weight k — 0.5). The jet charge is defined 
as: 

where p\ is the longitudinal momentum component with respect to the jet axis and qi is the charge of track 
i. 

2. Vertex charge, Q v t x = ^iWiQi, where u>i is the weight for track i to have come from a secondary instead of 
a primary vertex. Weights are derived from an artificial neural network with three inputs: the momentum 
of track i, the transverse momentum of track i with respect to the jet axis, and the impact parameter. 

3. Uncertainty on the vertex charge, crQ„ tx = T<iWi(l — Wi)q? . 

4. The product of lepton charge and output of neural network used to select b — > 1 decays, a lepton tag. 

The jet charge and vertex charge are not charge symmetric due to detector effects resulting in the difference in the 
rate and reconstruction of positively and negatively charged tracks. Offsets are evaluated using inclusive samples of 
T-tagged events and are subtracted prior to constructing Qt- The lepton tag is diluted by b mixing, cascade b decays, 
c decays, and fake leptons. Separate networks with fewer inputs are trained for events with a vertex or a lepton only. 

The variable Qt is defined as: 

_ N b (x) - N b (x) 

N b (x)+N b (x) [i) 
|Qt| = 1-2*7, (8) 

where Nb(x) and Ng(x) are the numbers of MC b and b hadron hemispheres with a given value of artificial neural 
network output x and r\ is the mistag probability. If Qt > 0, then the tagged hemisphere is more likely to contain a 
b hadron than a b hadron, and vice versa. If Qt = 0, both hypotheses are equally likely. 

The Qt flavor tag has some sensitivity to the decay flavor of the tagged hemisphere, which is not desirable in an 
inclusive measurement. Therefore, another tag Qm is applied in the opposite, or measurement, hemisphere. The Qt 
output in the T-tagged hemisphere, as well as the jet charge in the measurement hemisphere, Qm (with momentum 
weight k = 0), are combined to construct the composite variable: 
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Q-2 



(1-Qt)(1 + Qm) 



(1 - Q T )(1 + Qm) + (1 + Q T )(1 - Qm) 



- 1. 



(9) 



Again, if Q2 > (Q2 < 0), the so-called "M-tagged" hemisphere contains a b-hadron tag (b-hadron) tag. The Q2 
variable is designed to be sensitive to the production, but not the decay flavor of the b-hadron, thus avoiding biases 
to the reconstructed proper time measurement. After flavor tagging, 394119 events remain in the data sample. 

B. Proper Decay Time Reconstruction 



The CP-violating parameter a cp can be extracted from the rate asymmetry distribution, A(t), as defined in Equa- 
tion ^. This is accomplished by calculating the number of b-hadron M-tags minus the number of b-hadron M-tags 
in bins of reconstructed proper time t and performing a binned x 2 fit to measure a cp . The b-hadron proper time is 
defined as: 



t 



mi,L 



72 " 



(10) 



where L is the hadron decay length, Ej, is the b-hadron energy, and nib is the mass of the b-hadron, taken be that of 
the B+ and B° (5.279 GeV) @. 

The hadron decay length, L, is reconstructed in the measurement hemisphere by first forming a "seed" secondary 
vertex using the two tracks with the largest impact parameter, do, relative to the primary vertex in the highest energy 
jet. All tracks with p > 0.5 GeV, |g?o| < 1 cm, and < 0.1 cm, which are consistent with the "seed" vertex are then 
added to it via an interative procedure. The secondary vertex must contain at least 3 tracks and have an invariant 
mass exceeding 0.8 GeV, assuming all constituent tracks are pions. To further eliminate badly reconstructed or fake 
secondary vertices, the secondary vertex must be kinematically consistent with a long-lived particle originating from 
the primary vertex. Secondary vertices meeting the above criteria are identified in approximately 70% of M-tagged 
hemispheres for both signal and background. The decay length L between the primary and secondary vertices is then 
calculated using the jet axis as a constraint. 

The b-hadron energy is computed by first estimating the energy of the b-jet by treating the event as a 2-body 
decay of a Z° into a b-jet of mass nib and another object. The charged and neutral fragmentation energy, Ebfrag, 



was estimated using the procedure described in |13|, involving the charged track weights w;, and the unassociated 
electromagnetic calorimeter clusters weighted according to their angle with respect to the jet axis. The b-hadron 
energy is then, E b = E bjct - E b fra g - 



The reconstructed proper time distribution described by Equation 10 is convolved with 2 Gaussians to account for 
detector resolution effects. The RMS widths of the resolution functions are 0.33 and 1.3 ps and are determined from 
Monte Carlo studies. About 65% of events lie within the narrower Gaussian. These resolution functions represent an 
average over all true decay proper times t. The non-Gaussian effects apparent in small slices of t due to contamination 
from primary vertex tracks are not critical in this analysis, as the result is not particularly dependent on accurate 
decay time resolution. 



III. FIT TO RE(eb) 



For each of 34 time bins i (in the range -2 to 15 ps), the asymmetry is calculated in 10 bins j of IQ2I 

N b - - N l - 

lJ <\Q2\> ij (N! j +Nf j y 



(11) 



and the error cr^oba is given by: 
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where A 7 ^ (iVy) is the number of events with Q 2 > (Q2 < 0). The factor 1/ < |Q 2 | >y corrects for the tagging 
dilution (mis-tagging), which reduces the observed asymmetry for imperfectly tagged events. 

The 10 estimates of A?: bs in each bin i are then averaged, weighting by (<7 A obs) -2 to get A? bs . A binned x 2 -fit to the 
reconstructed proper time which accounts for non-B° background, the B° lifetime, and the Gaussian time resolution 
functions yields: 



a cp = 0.005 ± 0.055 ±0.013 
Re{e B ) = 0.001 ± 0.014 ± 0.003. 



(12) 
(13) 



The asymmetry, A(t), as a function for reconstructed proper time t is shown in Figure 0. The dots denote OPAL 
data, the solid line the fit result, and the dashed line the expected asymmetry for a cp = 0.15. 
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FIG. 1. Asymmetry of tagged b and b- hadrons as a function of reconstructed time t in data (points), and fit (solid line). 
The expected asymmetry for a cp = 0.15 is shown by the dotted line. 

The systematic uncertainties are summarized in Table ^. Detailed descriptions of the various contributions can be 
found in Reference |L4| . If the reconstruction efficiency for B° decays to different numbers of charm hadrons is not 
the same, the expected asymmetry could take the form: 



A(t) = c cp sin (Arridt) — a cp sin 2 (Arridt/2) 



Repeating the fit, letting both a cp and c cp vary gives: 



0.002 ±0.055 
0.026 ±0.027. 



(14) 



(15) 
(16) 



Differences in efficiency are not significant, as a cp does not change much. The systematic uncertainties on the 
measurement of c cp are listed in the second column of Table || 



IV. FIT TO (At/t) b 

The fractional difference between the b and b- hadron lifetimes is measured by dividing the data into 20 bins of Q2 
(related to b/b-hadron purity) and performing a simultaneous fit of the reconstructed proper time distrbutions to the 
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TABLE I. Systematic uncertainties on the measurements of a cp and c, 



Source 


Aa cp 


Ac cp 


B u lifetime 


0.002 


0.000 


Amj value 


0.001 


0.001 


h> traction 


0.002 


0.002 


Flavor tagging offsets 


0.003 


0.013 


Flavor tagging mis-tag 


0.009 


0.005 


b fragmentation 


0.008 


0.006 


Time resolution 


0.002 


0.000 


Total 


0.013 


0.015 



expected proper time distribution. The fit yields values for both T avg and (Ar/r)b- However, the result for T avg has 
a large systematic uncertainty due to the time resolution function and should not be interpreted as a measurement 
of the average b-hadron lifetime. The expected distribution accounts for time resolution effects, non-bb background, 
the lifetimes of b and b-hadrons, and a background component with a lifetime Tb g . The fit result is: 

(Ar/r) 6 = -0.001 ±0.012 ±0.008. (17) 

The uncertainty in the flavor mistag rate dominates the systematic uncertainty on (Ar/r)b. Reconstructed proper 
time distributions in 4 ranges of Q2 are shown in Figure 0. 



V. CONCLUSION 



An inclusive sample of b-hadron decays is used to search for CP and CPT violation effects. No such effects are 
seen. From the time dependent asymmetry of inclusive B° decays, the CP violation parameter is measured to be: 

Re{e B ) = 0.001 ±0.014 ±0.003. (18) 

This result agrees with the OPAL measurement using semileptonic b decays: Re(eB) = 0.002 ± 0.007 ± 0.003 0, and 
is also in agreement with other less precise results from CLEO and CDF || . The fractional difference in the lifetimes 
of b and b-hadrons is also measured to be: 

(Ar/r) b = -0.001 ±0.012 ±0.008. (19) 

This is the first analysis accepted for publication which tests the equality of the b and b-hadron lifetimes. These 
results are summarized in Figures I and |. 
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FIG. 2. Reconstructed proper time distributions in four ranges of Q2. The data are shown by the points with error bars and 
the fit prediction by the solid lines. The expected distributions for (Ar/r)b = 0.2 are shown by the dotted lines. 
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